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The apparent difficulty in recovering classical nonlinear dynamics and chaos from standard quan- 
tum mechanics has been the subject of a great deal of interest over the last twenty years. For open 
quantum systems - those coupled to a dissipative environment and/or a measurement device - it has 
been demonstrated that chaotic-like behaviour can be recovered in the appropriate classical limit. 
In this paper, we investigate the entanglement generated between two nonlinear oscillators, coupled 
to each other and to their environment. Entanglement - the inability to factorise coupled quantum 
systems into their constituent parts - is one of the defining features of quantum mechanics. Indeed, 
it underpins many of the recent developments in quantum technologies. Here we show that the en- 
tanglement characteristics of two 'classical' states (chaotic and periodic solutions) differ significantly 
in the classical limit. In particular, we show that significant levels of entanglement are preserved 
only in the chaotic-like solutions. 
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The correspondence principle is one of the fundamen- 
tal building blocks of modern physics. Stated simply, 
the principle insists that any new theory should contain 
the existing (experimentally verified) theories in some 
appropriate limit or approximation. Without this cor- 
respondence, it would not be possible to build new phys- 
ical theories without undermining previous theory and 
experiment. When quantum mechanics was originally 
proposed, the natural expression of the correspondence 
principle could be stated as '// a quantum system has a 
classical analogue, expectation values of operators behave, 
in the limit ft — > 0, like the corresponding classical quan- 
tities ' Whilst this works for many systems, problems 
sometimes arise when dealing with nonlinear systems . 
More generally, there are often found to be difficulties 
obtaining nonlinear classical equations from the linear 
Schrodinger equation. This is a significant problem and 
a great deal of work has been done to analyse the way in 
which classical nonlinearities express themselves in iso- 
lated (Hamiltonian) systems. In this paper we deal with 
an associated problem, namely nonlinear behaviour of 
quantum systems that are coupled to an environment 
(open quantum systems). 

Three main approaches exist for dealing with open 
quantum systems. The first, based on work by Feyn- 
man and Vernon treats the environment as a source 
of dissipation and decoherence and then calculates the 
average evolution of the quantum system by integrating 
out the environmental degrees of freedom. This provides 
a good method for large ensembles of quantum systems, 
but the average evolution does not reflect all types of non- 
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linear behaviour, such as chaos. The second approach is 
to examine the evolution of the probability distribution 
of the quantum system in phase space, having averaged 
out the environment Q. In this case, it can be shown 
that the quantum evolution of the probabilities is very 
similar to the evolution of the corresponding classical 
system, which is also coupled to a noisy environment. 
The third approach, often called quantum trajectories 
H HUSH is the one adopted here. This method 
uses continuous weak measurements on the environment 
and results in modified quantum evolution, which can 
be described by a stochastic Schrodinger equation. This 
equation reduces to the previous cases when averaged 
over an ensemble, but provides different trajectories for 
individual systems. The type of trajectories depend on 
the type of measurement that is performed on the en- 
vironment 0. Different trajectories are often referred 
to as unravellings of the Master equation, which would 
describe the average quantum evolution. For chaotic sys- 
tems, quantum trajectories show behaviour that is very 
similar to their classical counterparts in the correspon- 
dence limit ft — > 0, 0, ^3 ■ This paper demonstrates 
that the chaotic-like quantum dynamics have other prop- 
erties that do not reduce to the classical limit when ft 
becomes small. We show that the entanglement between 
two coupled chaotic oscillators can be significant - even 
when the trajectories are essentially classical. 

In this paper we have chosen to consider, as a con- 
venient example, two coupled Duffing oscillators that 
arc described by Quantum State Diffusion (QSD) 0, Q- 
Duffing oscillators were chosen because they are a stan- 
dard example of classical chaos, describing a driven oscil- 
lator that contains a third order nonlinearity. In addition 
single Duffing oscillators have been studied extensively 
using quantum trajectory models 0, 0|. The QSD 
unravelling corresponds to a unit-efficiency heterodyne 
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FIG. 1: The dynamics of the expectation values of position and momentum as a function of time (normalised to drive periods) 
for (3 — 0.01, /3 = 0.1, /3 — 0.25 and /3 = 1.0. The dynamics have been taken over the same duration and displayed so that, 
with the exception /? — 1.0, they approximately align at the time (re) at which entrainment occurs. 



measurement (or ambi-quadrature homodyne detection) 
on the environmental degrees of freedom 0. Similar 
chaotic-like dynamics can be obtained from the Duff- 
ing oscillator using most unravellings, and we consider 
another unravelling {quantum jumps - corresponding to 
photon number detection on the environment) below, and 
we show that similar entanglement properties can be ob- 
tained. In QSD, the evolution of the state_vector is 
given by the increment (Ito equation) \~ 
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where the Linblad operators Li represent coupling to the 
environmental degrees of freedom, dt is the time incre- 
ment, and are complex Weiner increments such that 
d|2 ^ ^ ^ and did^ = II 0- The first term 
on the right hand side of deals with the Schrodinger 
evolution of the system while the second (drift) and third 
(fluctuation) terms describe the decohering effects of the 
environment on the evolution of the state vector. To 
solve these equations can be a computationally demand- 
ing problem and approaching the classical limit generally 
requires the use of many basis states and extremely long 
run times. 



In the classical limit, QSD for one Dufflng oscillator 
reproduces the chaotic attractors (Poincare sections) of 
the classical equations 0, ■ this work, we extend 
this analysis and consider two identical, coupled Duff- 
ing oscillators. The Hamiltonian (total energy) for each 
oscillator is given by 
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where qi and pi are the position and momentum opera- 
tors for each oscillator. Here the Linblad operators are 
simply Li = yJ2Tiai (for i = 1,2), where is the os- 
cillator lowering operator. Following [T^ IT^ , we set the 
parameters for the drive and the damping to be gi = 0.3 
and Fj = 0.125 respectively. The Hamiltonian for the 
coupled system then takes the form 



H = Hi + H2 + ^.qiq2 



(3) 



where /i = 0.2 is a measure of the strength of the coupling 
between the two oscillators. 

Now, if we are to consider the correspondence prin- 
ciple in more detail we see that ft — !■ is not the only 
interpretation. By this we mean that, mathematically, 
it is equivalent to consider h remaining fixed (as indeed 
it does) and scale the Hamiltonian so that the relative 
motion of the expectation values of the observables be- 
comes large compared with the minimum area {h/2) in 
the phase space. In this paper, this is achieved by the 
introduction of a scaling parameter jS, where the system 
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behaves more classically as (3 tends to zero from its max- 
immn value of one. The choice of gi^Ti and the various 
values of (3 come from the previous work on the QSD 
Dufhng oscillator and dissipative quantum chaos |l2lll3l |. 

Wc solve the equation of motion for the state vector 
(|T|) using the Hamiltonian (j^J starting with the initial 
states centred on points taken from the Poincarc section 
of the uncoupled classical oscillators. Explicitly, this ini- 
tial state is a tensor product of coherent states for which 
the expectation values in position and momentum are 
centred in q — p phase plane at points chosen from the 
classical Poincare section. In Figure 1, we show the evo- 
lution of the expectation values of position and momen- 
tum for the two coupled oscillators (in blue and red) as 
a function of time for a selection of values of j3. From 
these figures it is apparent that the initial phase space 
trajectory of the coupled system is. for low (3 - the classi- 
cal limit, very similar to the that of the classical DufRng 
oscillator. In this small (3 limit, the uncertainties in posi- 
tion and momentum are a very small fraction (typically 
a few percent) of the oscillations, indicating that the dy- 
namics are essentially classical. As (3 increases to unity, 
this structure becomes harder to resolve and finally in 
the fully quantum limit (/? = 1) it is indiscernible. How- 
ever, it is clear, for (3 = 0.01,0.1 and 0.25, that after a 
certain time period the chaotic behaviour of the coupled 
oscillators gives way to regular, almost periodic motion. 
The oscillators can be said to have entrained^ i.e. their 
nonlinear oscillations are stabilised and synchronised by 
the coupling. We note for future reference that for (3=1, 
it is impossible to determine the presence of a chaotic or 
periodic attractor. Also, it is not possible to determine 
whether or not there is any level of entrainment between 
these oscillators. This is because, in the quantum limit, 
the quantum noise is large enough to drive the system 
away from periodic attractors and to swamp the fractal 
structure found in the chaotic classical attractor. 

Wc note that, when the two oscillators are entrained, 
the classical correlations between the oscillators are at a 
maximum (the oscillations arc, after all, fairly well syn- 
chronised). It is natural, therefore, to ask what level of 
quantum correlation is present within this system and 
how the chaotic and periodic phases of motion affect this 
correlation. We measure this quantum correlation by 
computing the entanglement entropy (mutual informa- 
tion) for the system from 

E = S{p,) + S{p2)~S{p) (4) 

where p = is the density operator of the system 

and Pi = Trj^i[/5] is the reduced density operator for 
each oscillator. Here S (pi) = Tr[pi In pi] is the von Neu- 
mann entropy |l3|. As QSD models the evolution of a 
pure state, S (p) = at all times. Since entanglement 
is a purely quantum mechanical phenomenon, one might 
expect that when a quantum system approaches the clas- 
sical limit the two oscillators would not be strongly en- 
tangled. 

In Figure 2 we show the entanglement entropy as a 
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FIG. 2: Entanglement entropy as a function of time (nor- 
malised to drive periods) for /3 = 0.01, /3 = 0.1, P — 0.25 and 
/3 — 1.0. By comparison with Figure 1 we can see that there 
is a sudden drop in entanglement entropy after the time Te at 
which the oscillators have entrained. 



function of time for the four different values of (3 as used 
in Figure 1, over the same corresponding time intervals. 
By comparing Figure 2 with Figure 1, we can see that 
its entanglement remains relatively large as long as the 
system undergoes chaotic motion. It is only when the 
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system becomes entrained that the entanglement drops 
approaching the classical limit. As we have already said, 
the uncertainties in the (3 = 0.01 case are very small 
when compared to the system's evolution in phase space. 
In this case, the system is well localised in phase space. 
The correspondence principle would seem to imply that 
these objects are behaving classically even though the 
entanglement (in the chaotic phase) does not seem to be 
diminishing. Once the oscillators entrain, the system be- 
comes approximately separable and the entanglement en- 
tropy falls as (3 reduces. These results would indicate that 
quantum oscillators operating in an apparently classical 
regime can generate significant levels of entanglement, 
but only in chaotic-like states. It is worth mentioning 
that the close connection between chaos in quantum sys- 
tems and entanglement has been noted by Wang et al. 
for Hamiltonian systems [Tsl l , and they have recently con- 
sidered the behaviour in the classical limit However, 
the system discussed in reference |l6j is non-chaotic and 
only entanglement between the system and its environ- 
ment is considered. Here, we consider the entanglement 
between two systems (each of which is coupled to an en- 
vironment). 
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FIG. 3: Mean entanglement entropy as a function of /3 for the 
chaotic-like and periodic (entrained) states. Here we see that 
the entanglement entropy for system in the chaotic state does 
not vanish as (5 approaches the classical regime. 

Each solution computed using quantum state diffusion 
corresponds to one realisation or experiment. We obtain 
a more general picture if we compute the average entan- 
glement as a function of (3 over many experiments: the 
mean entanglement entropy [l5l |. rather than the entan- 
glement of the average (i.e. mixed) state. We show the 
results of such a computation in Figure 3. The mean en- 
tanglement in the chaotic state does not appear to vanish 
as /? — > and, possibly more surprisingly, its maximum 
is not at /3 = 1. However, we also see that when this 
system entrains (at smaller (3 values) the entanglement 
is rapidly suppressed. Figure 3 also shows lines indicat- 
ing the best fit to the mean data points. In the case of 



the entrained state, the best fit is an exponential, which 
includes the quantum limit (/3 = 1), even though there is 
no evidence for an entrained state for this /3-value. 

For the preceding results, we have used the quan- 
tum state diffusion unravelling of the master equation. 
This unravelling makes assumptions about the underly- 
ing measurement process which is applied to the envi- 
ronment: unit efficiency heterodyne detection for QSD. 
However, studies have shown that the level of entangle- 
ment is dependent upon which unravelling is used |l7j . 
It is natural then to ask whether the properties of the en- 
tanglement that we have described above are particular 
to the QSD unravelling. To consider this, we choose an- 
other unravelling - the quantum jumps unravelling 0,0. 
This model is very different from QSD as it is based on 
a discontinuous photon counting measurement process, 
rather than a diffusive continuous evolution. However, 
it has already been shown that this rather different un- 
ravelling can produce, in the classical limit, the expected 
chaotic dynamics of the Duffing oscillator K seems 
natural, therefore, to use the quantum jumps unravelling 
as a second test unravelling for our example system. In 
this model, the pure state stochastic evolution equation 
for quantum jumps is given by 
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where dNj is a Poissonian noise process such that 
dNjdNk = SjkdNj, dNjdt = and dN~ = (l]Lj^ dt, 
i.e. jumps occur randomly at a rate that is determined 
hy{L]L,). 

As we can see from figure01the quantum jumps unrav- 
elling produces entanglement that persists in the classical 
limit so long as the oscillators remain in the chaotic-like 
solution. Once again, if the oscillators are seen to entrain, 
then the entanglement entropy drops significantly. The 
dependence of the entanglement entropy on (3 is similar 
in form for both unravellings (one continuous and one 
discontinuous). We find, as might be expected f^at 
the entanglement entropy for quantum jumps differs from 
that of QSD. However, even though the average level dif- 
fers between unravellings, the behaviour of the average 
entanglement - in both chaotic-like and non-chaotic - as a 
function oi (3 is very similar. We note that the data points 
shown in figure 0] exhibit a slightly more erratic nature 
than those shown in figure|31 The explanation for this lies 
in the fact that the the entanglement entropies obtained 
for the evolution of the system under the the quantum 
jumps unravelling take much more time to settle to a 
good average than for the QSD unravelling. Without 
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FIG. 4: The calculation of figure |S| using quantum jumps in- 
stead of quantum state diffusion. Again we show the mean 
entanglement entropy as a function of j3 for the chaotic-like 
and periodic (entrained) states. As with quantum state diffu- 
sion we see that when using quantum jumps the entanglement 
entropy for system in the chaotic state does not vanish as (3 
approaches the classical regime. 

very long computational times (e.g. several months) we 
cannot achieve much more significant accuracy with the 
computational power currently available to us. 



In this paper we have considered the classical limit of 
two coupled nonlinear quantum oscillators. The analysis 
is based on standard examples from classical chaos and 
open quantum systems. We have shown that when these 
oscillators enter the almost periodic entrained state the 
entanglement falls rapidly as the system approaches the 
classical regime, as expected. Conversely, we have also 
shown that these oscillators can remain significantly en- 
tangled even when the dynamics of the system appear 
to be classical, although only when in the chaotic state. 
Moreover, by using both the quantum state diffusion and 
quantum jumps models, we have shown that these results 
are not unique to a particular unravelling of the master 
equation. We find the results presented in this paper are 
surprising, since the quantum correlations represented by 
the non-zero entanglement do not have classical counter- 
parts yet they persist in the classical (correspondence) 
limit. 
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